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ABSTRACT 

We present the results of a ~230 ks long X-ray observation of the relativistic double-pulsar system 
PSRJ0737-3039 obtained with the XMM-Newton satellite in 2006 October. We confirm the detection in X- 
rays of pulsed emission from PSR J0737-3039A (PSR A), mostly ascribed to a soft non thermal power law 
component (F ^ 3.3) with a 0.2-3 keV luminosity of ^1.9 x 10^" erg s~^ (assuming a distance of 500 pc). 
For the first time, pulsed X-ray emission from PSR J0737-3039B (PSRB) is also detected in part of the orbit. 
This emission, consistent with thermal radiation with temperature kgT ~ 30 eV and a bolometric luminosity of 
'--^10^^ erg s~^, is likely powered by heating of PSRB's surface caused by PSR A's wind. A hotter (^130 eV) 
and fainter (^5 x 10^ erg s~') thermal component, probably originating from back-falling particles heating 
polar caps of either PSR A or PSR B is also required by the data. No signs of X-ray emission from a bow shock 
between PSR A's wind and the interstellar medium or PSRB's magnetosphere are present. The upper limit on 
the luminosity of such a shock component (^^10^^ erg s~') constrains the wind magnetization parameter um of 
PSR A to values greater than 1 . 
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The short-period, double-radio pulsar system 
PSRJ0737-3039 (Burgay et al. 2003; Lyne et al. 2004), 
besides being of paramount interest as a probe for theories 
of strong field gravity (Kramer et al. 2006), represents a 
unique laboratory for studies in several fields, ranging from 
the equation of state of super-dense matter to magneto- 
hydrodynamics. The system is observed nearly edge-on and 
consists of a fast, recycled radio pulsar (PSR J0737-3039A 
[PSR A], period P = 22.7 ms) orbiting its slower companion 
(PSRJ0737-3039B [PSRB], P = 2.77 s) with an orbital 
period of only 2.4 hr. Radio observations of the two pulsars 
permit one to derive a wealth of information which is not 
available for other neutron stars. These obviously include 
the neutron star masses and all the other geometrical and 
dynamical parameters of the system, but also information 
on the structure and physical properties of the two mag- 
netospheres. The double-pulsar is rich in observational 
phenomena, including a short radio eclipse of A by B 
and orbital modulation of the radio flux of B due to the 
influence of A (Lyne et al. 2004). The individual pulses 
from B show drifting features due to the impact of the 
low-frequency electromagnetic wave in the relativistic wind 
from A (McLaughlin et al. 2004b), while the eclipse of A 
is modulated at half the rotational period of B (McLaughlin 
et al. 2004c). High-energy observations are an important 
complement to these studies, in particular for what concerns 
the physics of the magnetospheric emission and dissipative 
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shocks in the close environment of the two neutron stars. 

X-rays could be pulsed magnetospheric or thermal emis- 
sion from pulsar A, as seen for several other recycled pulsars 
(Zavlin et al. 2002), but they could also originate in the col- 
liding winds of A and B (Lyutikov 2004). Because of the 
interaction of the relativistic wind of the fast spinning pul- 
sar A {E^^ = 5.8 X 10^^ erg s~') with the magnetosphere of 
its much less energetic companion (Ef^^ = L6 x 10^" erg s~'), 
the formation of a bow shock, likely emitting at high energies, 
is expected. The predicted fluxes are roughly comparable to 
those observed for magnetospheric and/or surface emission 
from the pulsars. A termination shock between the two pul- 
sars can probe the properties of a pulsar's relativistic wind at a 
smaller distance from the central engine than ever studied be- 
fore. In addition, detection of an orbital phase dependence in 
the X-ray emission might be expected (Lyutikov 2004; Arons 
& Tavani 1993). Such variability could constrain the geome- 
try of the emission site, thus providing new insights into the 
wind physics close to the pulsar Alternatively, most of the 
high-energy emission could arise from the shock generated 
when one or both the pulsar winds interact with the interstel- 
lar medium (Granot & Meszaros 2004). 

The first X-ray observation of PSR J0737-3039, a short (10 
ks) Chandra pointing (McLaughlin et al. 2004a) yielded only 
'^SO photons. The low X-ray luminosity (Lx = 2 x IO^^c/qj 
erg s~', where we indicate with d^ the distance in units of 
kpc) corresponds roughly to the entire spin-down luminosity 
of the slow pulsar B and to only a small fraction of that of 
pulsar A. The Chandra data of this pioneer observation could 
only poorly constrain the source spectrum, which appeared 
quite soft, and cannot provide significant evidence for vari- 
ability, due to the small statistics and time resolution being 
insufficient to look for the spin periods of the two pulsars. 

A longer public observation (50 ks), carried out in 2004 
March through the XMM-Newton Director's Discretionary 
Time program, yielded an improvement in statistics by a fac- 
tor ^10 with respect to the first look of Chandra. These data 
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confirmed the softness of the spectrum, which could be fit 
either by a power law with photon index F = 3.5;!;q 3 and ab- 
sorption ATh ~ 7.0 X 10^0 cm or by a blackbody with tem- 
perature of 0.15 keV and a lower interstellar absorption (Pel- 
lizzoni et al. 2004). These authors could also have performed 
the first X-ray timing analysis for PSR J0737-3039, but no 
periodic or aperiodic variations were found, with upper limits 
of ~60% on the pulsed fractions of both pulsars and of ~40% 
for modulations in the orbital period (all limits are at the 99% 
confidence level and for sinusoidal light curves). Campana 
et al. (2004) reported a joint spectral analysis of the XMM- 
Newton (only the MOS) and Chandra data, reaching similar 
conclusions on the source spectrum. 

A further Chandra observation ('^90 ks) with the High- 
Resolution Camera showed X-ray pulses at the period of 
PSR A (Chatterjee et al. 2007), with a double-peaked pro- 
file, similar to that observed in radio, and a pulsed fraction 
of ^70%. Although purely non thermal emission is consis- 
tent with the data, the X-ray pulse morphology of PSR A, in 
combination with previously reported spectral properties of 
the X-ray emission, suggests the existence of both non ther- 
mal magnetospheric emission and a broad sinusoidal thermal 
emission component from the neutron star surface. No pul- 
sations were detected from pulsar B, nor evidence for orbital 
modulation. Here we report the results of a ^230 ks long X- 
ray observation of PSR J0737-3039 obtained within the frame 
of XMM-Newton "Large Programs" in 2006 October. 

2. X-RAY OBSERVATION AND DATA REDUCTION 

Our observation of PSR J0737-3039 was carried out in two 
consecutive XMM-Newton orbits. The first part of the ob- 
servation started on 2006, October 26 at 00:28:15 UT and 
lasted 119.8 ks, the second one started on October 28 at 
00:28:37 UT and lasted 114.5 ks. In total, the observation 
allowed the coverage of ~26 revolutions of the binary system 
(Porb = 8, 834.535 s). The pn camera (Struder et al. 2001) was 
operated in Small Window mode (imaging across a 4' x 4' 
field of view with a 5.67 ms time resolution) with the medium 
optical filter. The MOS cameras (Turner et al. 2001) were 
also set in the Small Window mode (yielding a 2' x 2' field 
of view in the central CCD, with a time resolution of 0.3 s) 
with the medium filter. The data were processed using stan- 
dard pipeline tasks (emproc and epproc) of the XMM-Newton 
Science Analysis Software version 7.1.0. 

In view of the faintness of PSR J0737-3039 in soft X-rays, 
particular care in selecting source photons and reducing back- 
ground contamination is crucial. We used only photons with 
pattern 0-12 for the MOS cameras, while for the pn we used 
pattern 0^ for energies above 0.4 keV and pattern for the 
E <0.4 keV energy range. Such an event selection allows 
one to reduce by a factor ^3.5 the background count rate in 
the 0.15-0.4 keV energy range in the pn camera, while leav- 
ing almost unchanged the source count rate. 

We evaluated the optimal selection of source events by 
maximizing the source signal-to-noise ratio in the 0.15-10 
keV range, as a function of: 

Source extraction region. — We considered different extrac- 
tion radii in the 10"-30" range; the background was extracted 
in each case from two rectangular regions located at the same 
distance from the readout node as the source region. 

Threshold for high particle background screening. — The 
observations are affected by a few short soft proton flares. 
We extracted a Ught curve in the 0.15-10 keV range with a 



TABLE 1 

Photon Harvest from the XMM-Newton Observations of 
PSRJ0737-3039. 



Total Background Exposure 

Orbit Detector Counts" fraction'' (%) time'^ (ks) 



1260 pn 2483 24.5 80.7 

1260 MOSl 652 21.7 111.5 

1260 MOS 2 750 26.1 112.3 

1261 pn 2222 23.6 74.9 

1261 MOSl 587 19.7 102.9 

1261 MOS 2 697 29.5 103.9 



Note. — See text for details on data reduetion and background screening. 
^ Total number of counts inside the source extraction region. 
^ Contribution of background to total number of counts. 
Good observing time after dead-time correction and screening for soft-proton flares. 

10 s time bin from the whole field of view for each camera. 
Following the prescription by De Luca & Molendi (2004), 
we identified the quiescent average count rate and considered 
different thresholds in the 2-8 a range from this average rate. 

The optimal choice turned out to be a source extraction ra- 
dius of 18" for the pn and the MOS 1 cameras and of 15" 
for the MOS 2 camera (which has a higher low-energy back- 
ground) and a threshold at 5(t from the quiescent rate to 
screen from high particle background episodes. The resulting 
source and background count statistics are shown in Table 1. 

3. TIMING ANALYSIS 

3.1. Pulsations from PSR J0737-3039A 

For the search of pulsations at the spin period of pulsar 
A we could only use the pn data (~4700 counts), owing to 
the inadequate time resolution of the MOS data. The times 
of arrival were converted to the solar system barycenter and 
corrected for the orbital motion and relativistic effects of the 
binary system according to Blandford & Teukolsky (1976). 
This was done with a program we wrote ad hoc and tested 
using an XMM-Newton observation of the binary millisec- 
ond pulsar XTE J1751-305 (Miller et al. 2003). As a further 
cross-check, we also compared the rotational phases obtained 
by our program with those resulting from the Tempo timing 
analysis software.^ The maximum discrepancy between tim- 
ing corrections from our program and Tempo is at most a few 
microseconds. 

In our search for pulsations we employed the test (Buc- 
cheri et al. 1983), where we indicate with n the number 
of harmonics. We examined a wide frequency range cen- 
tered on the value predicted at the epoch (MJD 54034) of 
our XMM-Newton observation by the radio measurements of 
Kramer et al. (2006). The most significant statistics oc- 
curred for n = 1 at P^^^ = 22.6993787(5) ms (1 a errors in 
the last digit are quoted in parentheses). The correspond- 
ing Zj value is 378.91, which even when taking into account 
the ^10"^ searched periods, has a virtually null probability of 
chance occurrence. Our best-fit period is consistent with the 
value = 22.6993784661 12(5) ms expected from the ra- 
dio ephemeris of Kramer et al. (2006). 

The background-subtracted and exposure-corrected light 
curve of pulsar A in the 0. 15^ keV energy range is shown in 
the top panel of Figure 1 (due to the soft pulsar spectrum, the 
signal-to-noise ratio above 4 keV is very small). Error bars are 

' See http://www.atnf.csiro.au/research/pulsar/tempo. 
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Fig. 1. — Background-subtracted pn light curve of PSR A in the 0.15- 
4 keV energy range folded at period P^'" = 22.6993787(5) ms (top) and 
corresponding hardness ratio curve based on the soft (0.15-0.3 keV) and hard 
(0.8-3 keV) energy ranges (bottom). The dashed lines indicate the three phase 
intervals (A, B, and C) used for the phase-resolved spectroscopy. 

calculated according to the expression a, = {Ci+Bif'^yi^ /Ei, 
where Q, Bi and £, are respectively the total counts, back- 
ground counts, and exposure in the bin and / is the ra- 
tio between the source and background extraction area (/ ~ 
0.11). We checked that the pulsations are significantly de- 
tected down to the lowest energy bins covered by the pn in- 
strument; in fact, the Pearson statistics for a 10 bin light curve 
in the 0.15-0.2 keV range gives a reduced (hereafter x^) 
of 3.85, corresponding to a Aa detection. 

The pulse profile is double peaked and deeply modulated, 
with a pulsed fraction of 75.7% ±5.4% in the energy range 
0.15^ keV. To calculate the pulsed flux, we considered 
all the counts above the minimum of the light curve, us- 
ing the expression PF = Qot - nA'min and its associated er- 
ror (TPF = {Ctoi + rp-ali^yi'^ ~ n{Nj^^yi^, where Qot are 
the total counts, n is the number of bins in the Ught curve 
and are the counts of the minimum. The expression 
(Cmax-Cmin)/(Cmax + Cmm) uscd by Chattcrjce et al. (2007) 
does not account for the pulsed flux associated with secondary 
peaks; nevertheless their pulsed fraction value is in agreement 
with our result within \a. Both methods are "bin depen- 
dent," but reasonable different choices of the number of bins 
(i.e. n > 10) do not significantly affect the results. Note that 
the pulsed fraction upper Umit of 60% obtained in the short 



2004 XMM-Newton observation referred to a sinusoidal pro- 
file (Pellizzoni et al. 2004), and therefore it is not inconsistent 
with the present result. 

A simple hardness ratio analysis, based on the soft (5: 0.15- 
0.3 keV) and hard (H: 0.8-3 keV) energy ranges,^ indicates 
a softer spectrum in correspondence to the minimum in the 
folded light curve (see bottom panel of Figure 1). To inves- 
tigate the energy dependence of the pulsed fraction in more 
detail and as a cross-check of complementary phase-resolved 
spectral studies discussed below, we produced the folded light 
curves for four different energy ranges (0.15-0.3 keV, 0.3-0.5 
keV, 0.5-0.8 keV, and 0.8-3 keV) providing 800-900 source 
counts each (Figure 2). The corresponding pulsed fractions, 
plotted in Figure 3, increase from 50% to 90% from lower to 
higher energies. These values have a probability smaller than 
1 % of deriving from an energy-independent distribution. 

We searched for possible flux and pulse profile variations 
due to the pulsars' mutual interaction as a function of the or- 
bital phase. For example, if there is a bow shock structure 
as that described in Lyutikov (2005), associated with an un- 
pulsed variable component (see § 5), we would expect vari- 
ations in the total flux correlated with changes in the pulsed 
fraction. In Figure 4 (top) we report the pulsed fraction and 
total flux obtained from the folded 0.15-3 keV light curves 
integrated in four orbital phase intervals. The pulsed fraction 
varies in the 50%-80% range with the minimum correspond- 
ing to the superior conjunction of A (when A is occulted by 
B; dotted line). However, the values are consistent at the 5.5% 
level with a uniform distribution and there is no correspond- 
ing change in the total flux. A similar analysis dividing the 
orbit in eight parts (Figure 4, bottom) gave only marginal evi- 
dence of a ^15% variability (null hypothesis probability nuU 
hypothesis probability = 1.5%) in the pn flux. 

To calculate upper limits on orbital flux modulation, we 
considered the fraction of the counts above the minimum of 
the light curve with error evaluation similar to those described 
in this subsection. For time scales of ^15 min (10 bins) the 
1 a upper limits on variability (pn-nMOS data) is of 1 1 .5% for 
the 0.15-4 keV range and <20%-30% for the selected bands 
mentioned above. Longer time-scales of 0.5-1 hr imply upper 
limits <10% on all energy selections. The search for orbital 
and aperiodic variability, even when selecting the time inter- 
vals corresponding to the minimum in the PSR A folded light 
curve, does not improve the above upper limits. 

3.2. Pulsations from PSR J07 37-30398 

We searched for pulsations from PSRB using the same 
method and radio ephemeris reference as for PSR A (Kramer 
et al. 2006), but in this case, owing to the longer pulse pe- 
riod, we could also use the lower time resolution MOS data. 
The expected pulsar period at the epoch of our observation 
is /»^dio ^ 2.7734607024(7) s, wifli an uncertainty orders of 
magnitude smaller than the intrinsic pulse search resolution 
of our data set, \pI/Tob^ ~ 1.4 x 10"^ s. Thus, if the rota- 
tional parameters of the neutron star are stable (as suggested 
by radio ephemeris), a single-trial pulse search at Pq is ap- 
propriate. However, since we cannot in principle exclude that 
(unlikely) gUtches and/or significant timing noise occurred in 
the time between the radio and X-rays observations, we also 
performed a search in a range of periods around P^^^°. In 



* The hardness ratio is defined as {H — S)/ {H + S), where for each phase 
bin H and S are the background-subtracted count rates in the hard and soft 
ranges. 
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Fig. 2. — Background-subtracted light curves of PSRA for four energy intervals providing 800-900 counts each. The count rate of the minimum of the light 
curve is compatible with at high energies, while at low energies the presence of significant unpulsed flux is apparent. 



A 



A 

4 



□ 



11. 



0.2 L 
0.1 



A 



Energy (keV) 



Fig. 3. — Pulsed fraction of PSRA (crosses) as a function of energy (en- 
ergy intervals as in Figure 2). The plot also shows the predicted pulsed flux 
(triangles, stars, and squares) for the spectral models discussed in §§ 4 and 
5. 



any case, no significant pulsations were detected in the whole 
energy range nor in our canonical bands (0.15-0.3 keV, 0.3- 
0.5 keV, 0.5-0.8 keV, and 0.8-3 keV) providing 1200-1400 
source counts each. The upper Umit on the pulsed flux frac- 
tion is of 40%. 

Since the radio flux of PSRB is strongly modulated as a 
function of the orbital phase and is nearly disappearing for 
most of the orbit (Lyne et al. 2004; Burgay et al. 2005), 
we also performed an orbital phase-dependent pulse search. 



We folded the counts at pi^'^^° for four orbital phase inter- 
vals around the pulsars' conjunctions (A0 = [0.16-0.41] and 
[0.66-0.91]) and quadratures (A0 = [0.41-0.66] and [0.91- 
1.16]). As a reference, consider that orbital phase 0.28 corre- 
sponds to the superior conjunction of A. A positive detection 
of the PSR B pulsations in the pn data set is obtained for the 
phase interval corresponding to the ascending quadrature of B 
(A(/) = [0.41-0.66]). The Pearson statistics for the pulse his- 
togram with 20 bins (^250 mean source counts per bin) in the 
0.15-3 keV band provides a Xr of 2.4 (19 degrees of freedom 
[dof]) corresponding to a probability of only 0.06% that the 
profile is drawn from a uniform distribution (3.4(t). A fur- 
ther analysis with the test supports the detection; with the 
number of harmonics n being varied from 1 to 5, we found a 
statistically significant signal for n = A (Z^ = 34.55, 4.2 a) and 
5 (Zj = 35.48, 3.8cr). No significant detections are obtained 

for the other orbital phase intervals (Xr < 1-4, 19 dof). 

We tried to better constrain the orbital phase interval corre- 
sponding to a positive detection of PSRB by shifting and/or 
shortening the orbital phase intervals of the period search. 
However we did not find more significant detections, and 
considering the low- statistics pulsed data, we cannot derive 
tighter constraints. 

As mentioned above, the stability of the timing parameters 
of this neutron star is likely, but not guaranteed. Thus, in prin- 
ciple, a strict single-trial search at the predicted value p^dio 
might not be the optimal strategy. Furthermore, a blind search 
in period is an important cross-check to claim a robust pulsed 
detection. We therefore also searched for other peaks in the Z^ 
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Fig. 4. — Pulsed fraction and total flux of PSR A as a function of orbital 
phase integrated on ~40 min (top) and ~20 min {bottom) bins. The dotted 
vertical line corresponds to the superior conjunction of A (when A is occulted 
byB). 
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Fig. 5 . — Distribution of the Z4 statistics for the pn data of PSR J0737-3039 
during the orbital phase interval A</) = [0.41-0.66] (see § 3.2 for details). The 
vertical dashed line indicates the PSRB period expected from adopting the 
radio ephemeris of Kramer et al. (2006). 

distributions around We did not find more significant 

detections apart from values close (within the pulse search 
resolution) to f '^^lo (pjguj-g 5^ ^he best-fit pulse period is 
P^^'t = 2.773459(4) s providing a probability of only r^\(r^ 
(Z| = 42.63, 4.9 cr) that we are sampUng a uniform distribu- 
tion. 

The resulting pn light curves for the four phase intervals 
considered are shown in Figure 6. The pulse profile appears 



highly structured and with a pulsed flux fraction of 35 ± 15% 
in the energy range 0.15-3 keV. 

The detection of PSRB in MOS data is only marginal, and 
the analysis of merged MOS and pn data does not provide bet- 
ter detection significance (even when selecting counts which 
are in the lowest emission bins of the pulse profile of PSR A), 
probably because the 0.3 s time resolution of the MOS dilutes 
the sharp peaks of the light curve. Furthermore, the poor ef- 
fective area of the MOS at low energies (i.e. negligible in 
the 0.15-0.3 keV band) is not efficient for pulse searches in 
sources with very soft spectra. In fact, although PSR B is also 
marginally detected at £ > 0.5 keV (x^ = 2.1, 19 dof), most 
of the flux contributions to the main peaks of its light curve 
come from the low-energy band. 

In the hypothesis that PSRB is emitting X-rays only in the 
orbital phases in which we can firmly claim pulse detection, 
an orbital flux variability with magnitude comparable with 
the pulsed flux fraction should be present. In fact, there is 
marginal evidence that the pn total flux peaks around (j) ~ 0.5 
when B's pulses appear (Figure 4). It is worth noting that 
the sum of the pulsed emission of the two pulsars roughly ac- 
counts for the whole X-ray flux, not allowing for a significant 
amount of emission from other sources (e.g. shock emission 
from the pulsars' mutual interactions). 

4. SPECTRAL ANALYSIS 

Source and background spectra were extracted from the pn 
and MOS data using the regions and soft protons filtering cri- 
teria described in § 2. For what concerns the pn pattern se- 
lection, we used only single events (pattern 0) for the whole 
energy range in order to reduce the background in the low- 
est energy channels. We checked that fully consistent results 
were obtained including double events (pattern 0-4). 

Each spectrum was fitted with the XSPEC version 11.3.2 
software package using four models, aU modified by pho- 
toelectric absorption (phabs model with cross section from 
Balucinska-Church & McCammon 1992 and abundances 
from Anders & Grevesse 1989): power law (PL), black- 
body (BB), two blackbodies (BBc+BBh), and power law plus 
blackbody (PL+BB). The best-fit parameters for the pn spec- 
trum and for a simultaneous fit of the pn, MOS 1, and MOS 2 
spectra are reported in Table 2. 

In contrast with the previous short Chandra and XMM- 
Newton observations (McLaughlin et al. 2004a; Pellizzoni 
et al. 2004; Campana et al. 2004), single-component mod- 
els fail to properly fit the data. A blackbody model is com- 
pletely ruled out and only a very steep power law provides 
a barely acceptable fit (null hypothesis probability = 4%). 
Either a double-blackbody or a power law plus blackbody 
model are instead an adequate fit (null hypothesis probabihty 
> 50%). In the two-component models, the derived colunm 
density is also in better agreement with the value estimated 
from the radio pulsars' dispersion measure 

(A^H = 1.5 X lO^o 

cm~^, assuming 1 H for each electron; Campana et al. 2004). 
The power law photon index in the PL+BB model is stiU 
very soft. A hard power law with photon index F ~ 2, as 
expected for a shock emission, does not provide an accept- 
able fit. The observed flux in the 0.2-3 keV energy range is 
4.0 X 10"'^ erg cm~^ s~\ corresponding to an unabsorbed lu- 
minosity Lx = 2.2 X lO^^dl^ erg s"' for the PL+BB model. 
The radio dispersion measure distance of ~500 pc is based on 
a model for the interstellar electron density (Cordes & Lazio 
2002), and the parallactic distance ranges from 0.2 to 1 kpc 
(Kramer et al. 2006), implying uncertainties of a factor -^4 in 
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Fig. 6. — Light curves of PSRB folded at period F^^^^ = 2.773459(4) s as a function of the orbital phase. The pulsar is detected at the 4.9 ct level around the 
ascending node (orbital phase interval 0.41-0.66). 



TABLE 2 

Summary of the Spectral Results in the 0. 15-10 keV Energy Range. 



Instrument 


Model 


Nh 
(lO^^cm-^) 


r 


PL norm." 


(eV) 


Rbbh* 
(m) 


ksTEBc 
(eV) 


RbBc" 
(km) 


Xr (dof) 


pn 

pn+MOS 


PL 
BB 

BBh+BBc 
PL+BB 
PL 
BB 
BBh+BBc 
PL+BB 


*-"-0.7 



<0.9 
3,2+1-3 

^•"-0.6 


<0.5 

3.0±1.1 


3.4±0.1 

o O+0.3 
J-J_0.4 

3.4±0.1 

•3 1+0.1 

-'• -0.3 


9.4±0.6 

6.6! j;^ 
9.6±0.4 

6.5±1.0 


142 

280!5'' 
160l 

147 

290!^!] 
150±20 


149 

23+'2 

DU_3Q 

138 

201 
80±30 


114!^ 


"•''-'-0.03 
(1 21+004 


1.38 (50) 
3.29 (50) 
1.14(48) 
1.20(48) 
1.16(98) 
3.05 (98) 
1.00 (96) 
0.92 (96) 



Note. — Errors are given at the 90% confidence level and are not reported for the single-blackbody model because the tit is largely unacceptable. 
" In units of 10^ photons cm"^ s"' keV"', at 1 keV. 



Assuming a distance of 500 pc. 



luminosity estimates. 

The good timing resolution coupled to the simultaneous 
spectral capability of the pn data, together with the rela- 
tively large number of detected source photons, allowed us 
to perform for the first time a phase-resolved spectral analy- 
sis for PSR A. Three spectra were extracted by selecting the 
pn events (with pattern equal to 0) in the three phase intervals 
shown in Figure 1, corresponding to the off -pulse (phase A), 
the interpulse (phase B), and the main pulsation peak (phase 
C). The results of the fit of these three spectra with the four 
models used for the phase-averaged spectrum are shown in 
Table 3. A single-blackbody model does not provide an 
acceptable fit for any of the three spectra. We also note that, 



while the best-fit parameters of the main peak (phase C) and 
interpulse (phase B) are very similar to those of the phase- 
averaged spectrum, the fits of the off-pulse spectrum (phase 
A) with the two-component models converge to rather differ- 
ent values for some parameters. In particular, in the double- 
blackbody model, one of the two blackbodies has a tempera- 
ture ^10 times smaller (and a correspondingly larger emitting 
area) than in the other two phase intervals. Moreover, in the 
PLh-BB model the power law component is steeper, and in 
both cases, the column density is significantly larger. 

To better check if the source spectrum is actually changing 
with the phase of PSR A, we simultaneously fitted the three 
phase-resolved spectra with the same blackbody plus power 
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TABLE 3 

Results of the Phase-resolved Spectroscopy Performed with the pn Single Events in the 0. 15-10 keV Energy Range. 



Model 


Phase 




r 


PL norm.^ 




Rbbh'' 


*:b7bbc. 


KbBc'' 


X? (dof) 




interval 


(lO^^cm^^) 






(eV) 


(m) 


(eV) 


(km) 




PT 

r^ij 


A 
/\ 

B 


7 1+2.2 

4.9t\i 


T c+0.5 
J.J_() 3 

o O+0.3 


2 9+"-'' 










0.99 (23) 




c 


5.o;|i 


3.4±0.2 












1.11 (42) 


BB 


A 
B 
C 









118 
141 
146 


149 
138 
167 






2.57 (9) 
2.16(23) 
2.67 (42) 


BBh+BBc 


A 


Q n+3.2 
y.u_3(| 






140 ±20 


130±40 


29+:*2 


35 g+93.3 


0.94 (7) 




B 


< 1.9 








ll+'3 


12050 


n 1 o+^'o^ 

^'■i'^-0.02 


0.86 (21) 




C 


<2.3 






260±50 


-'^-12 


10050 


A ^n+0-31 


1.00(40) 


PL-l-BB 


A 


10±3 


q q+14.3 


0.02!»;™ 


140 ±20 


1 on+eo 






1.14(7) 




B 


^■"-1.8 


2.7 ±0.7 


5 1+2.8 


120ii« 


0.14±0.09 






0.88 (21) 




C 


4.0±1.7 


-'■-'-0.3 


^■'-3.2 


190!i 


50+50 






0.95 (40) 



Note. — Eirors are given at the 90% confidence level and are not reported for the single-blackbody model because the fit is largely unacceptable. 
^ In units of 10"*^ photons cm~" s~' keV~\ at I keV. 
^ Assuming a distance of 500 pc. 
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Energy (keV) 

Fig. 7. — Simultaneou.s fit of the three pn phase-resolved spectra (see Fig- 
ure 1 for the coiTesponding phase intervals) with a power law plus blackbody 
spectrum with linked parameters, except for an overall normalization factor 
(see Table 3). For the spectrum of the pulse peak (black), the relative con- 
tribution of the power law {blue) and blackbody {light blue) components are 
also shown in the spectrum. The residuals are shown separately for each 
spectrum in standard deviation units. 

law model with all the parameters linked together, except for 
an overall normalization factor A globally acceptable fit is 
obtained, but if we compute separately the Xr for each, fixing 
all the parameters to the best-fit values of the simultaneous 
fit, we note that the fit of the off-pulse spectrum is unaccept- 
able (see Table 4). The corresponding residuals (see Figure 7) 
show a soft excess in the off-pulse spectrum, confirming the 
softening in the pulse minimum already displayed in the hard- 
ness ratio (see Figure 1). 
We also tested the possibility that separate physical pro- 



cesses, producing different spectral components, were respon- 
sible for the pulsed and unpulsed emission. A first hypothesis, 
already proposed by Chatterjee et al. (2007), where a black- 
body accounts for the unpulsed emission and a power law for 
the pulsation is immediately discarded by the fact that the off- 
pulse spectrum cannot be well fitted by a single-blackbody 
model (see Table 3). 

The off-pulse spectrum is instead well fitted by two- 
blackbodies. Therefore we tried to fit the phase-resolved spec- 
tra with a three-component model composed of two constant 
blackbodies and a power law with stable photon index and 
variable normalization. The resulting fit is very good and its 
parameters are shown in Table 4. 

Since a good fit of the off-pulse spectrum is also obtained 
using a PLh-BB model (see Table 3), we also tested a model 
where the unpulsed spectrum was composed by a blackbody 
plus a power law and a second blackbody was responsible 
for the pulsed emission. In this model, the two blackbodies 
might come from the surface of the two pulsars and the power 
law from a shock. However, the resulting fit, also reported in 
Table 4, is not acceptable. 

5. DISCUSSION 

Our deep XMM-Newton observation of the double-pulsar 
PSR J0737-3039 reveals a complex X-ray phenomenology 
that cannot be simply ascribed to PSR A alone, although it 
is clear that the only significant power plant of the system 
(apart from its strong gravitational potential) is the spin-down 
energy release of this pulsar 

We confirmed deeply modulated pulsed X-ray emission 
from PSR A and we also detected for the first time X-ray pul- 
sations from PSRB. The rotational energy loss of PSRB is 
by far insufficient to produce the observed X-ray luminosity. 
It can also be excluded that the observed PSRB's X-rays de- 
rive from residual internal heat of this 50 million years old 
pulsar. Therefore, this detection gives strong evidence for the 
mutual interactions of the two pulsars also at high energies. 
Although, as discussed below, there is no direct evidence of a 
termination shock between their winds and magnetospheres, 
we interpret the PSRB emission, visible only at the orbital 
phases in which the PSR A beam intercepts PSRB, as ulti- 
mately powered by PSR A. 

We found that double-component models are required to 
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TABLE 4 

Results of Simultaneous Fit of the Three Phase-resolved pn Spectra of PSR A. 



Fixed 
component 


Variable 
component 


Phase 
interval 


(lO^^cm'^) 


r 


PL 
norm.^ 


feV") 


SbBh'' 

Cm) 


(eV) 


^BBc'' 

(km) 


-^r,single 

(dof) 


Xr,tot 

(dof) 


BBh+BBc 
PL+BBc 


PL+BB= 
PL 
BBh 


A 
B 
C 
A 
B 
C 
A 
B 
C 


3.2t[i 
<1.6 


.^•^-0.3 

■J o+O.l 
-'•-'-0.2 

O+0.5 
^•^-0.7 


^■"-0.8 

5.9±1.5 
9.8!J| 
<2.2 

D.U_i 4 

1 O + l.O 

^■•^-L3 


166±14 


45+30 

110±40 

<45 

70+20 


32!^ 


15!^ 
0.26!^;^5 


2.15 (12) 
0.84 (26) 
0.89 (45) 
0.71 (12) 
0.96 (26) 
0.91 (45) 
1.70(12) 
0.94 (26) 
1.25 (45) 


1.15 (76) 
1.00 (74) 
1.37 (74) 



Note. — Errors are given at the 90% confidence level. 
^ In units of 10"^' piiotons cm"" s"' IceV"^ at 1 keV. 
^ Assuming a distance of 500 pc. 

^ The relative normalization between the two components is fixed. 



account for the spectral phenomenology of the double-pulsar 
X-ray emission. A double-blackbody (BBc, ksT = 114 eV 
and BBh, ksT = 290 eV) provides a good fit, while a PLh-BB 
model with a blackbody temperature of fegT = 150 eV and a 
soft photon index of 3.2 is even better. 

It is interesting to compare the relative contribution of each 
spectral component to the total flux with the pulsed fraction 
of PSR A in order to test whether the pulsed flux could be as- 
cribed to a specific thermal or non thermal component (or to 
a particular combination of them). For the case of the PLh-BB 
model, the power law component contributes to ~76% of the 
pn counts in the 0.15^ keV band, a value in perfect agree- 
ment with the pulsed fraction of PSR A (75.7 ±5.4)% in the 
same energy range. However, the tentative association of 
most of the non thermal flux in the PLh-BB model with the 
pulsed flux is ruled out by observing that, in this hypothesis, 
the steep power law should strongly contribute to the pulsed 
flux especially at low energies, but this is not the case; the 
power law contributes ~86% of the counts in the 0.15-0.3 
keV band, a value much higher than the corresponding pulsed 
flux (~50%) as shown by the stars plotted in Figure 3. A 
similar argument applies for the case of the double-blackbody 
model (see triangles plotted in Figure 3). These results are 
confirmed by the phase-resolved spectra of PSR A showing 
that the off -pulse spectrum cannot be well fitted by a single- 
component model. 

More generally, none of the double-component best-fit 
models are compatible with the plot in Figure 3, even with 
the assumption that the pulsed flux results from a combina- 
tion of both spectral components. In fact, the total flux of the 
softer spectral component (either the power law or the colder 
blackbody) is 3^ times higher than that of the harder com- 
ponent. Therefore, it is impossible to invert the trend in Fig- 
ure 3 even when assuming that the harder component is fully 
pulsed. For example, in the case of the PL+BB model, as- 
suming that the pulsed flux is due to the full blackbody com- 
ponent plus a part of the power law component, it is impossi- 
ble to obtain a pulsed fraction that is increasing with energy. 
Therefore, the entire X-ray flux of the double-pulsar system 
cannot be interpreted as a simple combination of a non ther- 
mal and single-temperature thermal emission from PSR A or 
PSRB. Therefore, on the basis of these considerations and of 
the spectral results of § 4 we are led to consider models based 
on three spectral components. 

The inclusion of an additional cooler blackbody BBc (pos- 



sibly related to PSRB) in the PL+BB model can account for 
the small pulsed fraction at low energy. Three-component 
models where blackbody emission (BBh) is responsible for 
most of the pulsed emission from PSR A do not provide an 
acceptable spectral fit (see Table 4) and fail to match pulsed 
flux fractions in Figure 3. Assuming instead that most of 
the pulsed PSR A flux is due to the power law, it is possi- 
ble to reproduce the trend seen in Figure 3 (squares). This 
is also the best three-component scenario compatible with the 
phase-resolved spectrum of PSR A (see § 4, Table 4). It is 
worth noting that the phase-resolved spectrum of PSR A does 
not reveal any variation apart from a significant softening in 
the off-pulse phases, in agreement with the proposed three- 
component model. 

The best three-component model fitting phase-resolved 
spectrum of PSR A and pulsed fluxes of both pulsars is thus 
composed of a power law with photon index F ~ 3.3 and 0.2- 
3 keV luminosity LpL = 1.9 x 10^^ ergs"', a cooler (~30 
eV) blackbody with bolometric luminosity Lbbc = 3 x lO^' 
erg s~', and a hotter blackbody (~130 eV) with luminos- 
ity Z'BBh = 5.2 X 10^^ erg s~^ The emission from PSR A is 
mostly non thermal, its luminosity being equal to LpL plus 
the possible additional smaller contribution of LbBh^ ^a ^ 
2 X 10^" = 3 X IQ-'^Ea in die 0.2-3 keV energy range. 

Because of the softness of the spectrum, the efficiency of 
the conversion of rotational energy loss into X-ray luminos- 
ity is strongly sensitive to the energy range considered. The 
PSR A luminosity in the 0. 1-10 keV range matches ~ 10~^iiA, 
the value typically found for recycled radio pulsars that are 
detected at X-ray energies (Becker & Triimper 1999). As ex- 
pected from a mostly magnetospheric emission, pulse peaks 
are still relatively steep and far from thermal broad sinusoidal 
pulsations seen in other recycled pulsars with comparable spin 
parameters (Zavhn et al. 2002). Non-thermal deeply modu- 
lated components are instead seen in the pulsed flux of the 
most energetic recycled pulsars (e.g. PSRB1821-24). The 
peculiar nature of PSR A is also stressed by the fact that no 
other recycled pulsar with such a steep non thermal emission 
component is known. 

The observed equivalent emitting radius of BBh is one or- 
der of magnitude smaller than the polar cap radius of PSR A 
(~ 1 km), but modeling with a non-uniform temperature of the 
heated region (Zavlin et al. 2002) or invoking a partially filled 
polar cap (Harding & Mushmov 2002) could solve this dis- 
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crepancy. Thus, we can interpret BBh as due to the heating 
of A's polar cap from back-flowing particles in its magneto- 
sphere. The light curve of PSR A shows a bridge between the 
two peaks 2-3 times higher than the off-pulse flux, leading 
to the idea that both non thermal pulses emerge from a wide, 
center-filled thermal emission cone from the magnetic pole 
possibly associated to BBh- 

5.1. The nature of PSR B X-ray emission 

We detected PSR B during part of the orbit, around the as- 
cending node. Its 0.2-3 keV flux estimated from the pulsed 
light curve is ~1.2 x 10"'''^ erg cm~-^ s~^ and corresponds to 
(35 ± 15)% of the total flux in the orbital phase interval 0.44- 
0.61. The corresponding luminosity of 3.6 x 10^^ erg s~^ rep- 
resents ^20% of the spin-down energy loss, a value much 
higher than observed in all the other normal pulsars (Possenti 
et al. 2002). 

The PSR B emission is likely associated to the cooler black- 
body emission BBc (ksT ~ 30 eV), with the possible addi- 
tion of the weaker and hotter component BBh (if the latter 
is not associated to PSR A). In this case, the bolometric lu- 
minosity of PSRB can be calculated as Z-b — 4 x Lbbq — 
1.2 X 10^^ erg s~' assuming it is emitted during one-fourth 
of the orbit (Lbbc ^bbh ^re orbital phase-averaged lumi- 
nosities). This luminosity is surprisingly much higher than the 
X-ray emission from PSR A (assuming La — LpL = 1 .9 x 10^*^ 
ergs"'). Obviously, X-ray emission from PSRB can only 
be powered by an external source, i.e. the spin-down energy 
from PSR A. In this hypothesis, ^2% of the rotational energy 
loss of PSR A is converted into thermal radiation from PSR B. 

Zhang & Loeb (2004) explain the strong brightening of 
pulsed radio flux from PSRB in two portions of the orbit as 
episodes during which pairs from PSR A's wind flow into the 
open field line region of pulsar B and emit curvature radia- 
tion at radio frequencies within an altitude of ~10^ cm. Once 
pairs from A's wind leak into B's magnetosphere, they will 
stream all the way down, heating the surface of B, giving 
rise to thermal radiation. According to this model, the pre- 
dicted X-ray luminosity is Lx,b ~ 2 x 10^'r? (Ar2w,A/47r)~' 
erg s~^ , and the typical temperature of the polar cap of PSR B 
is ksTpc ^ 0.4r/i/4(Af7w,A/47r)-i/'^ keV, where rj is the frac- 
tion of all pairs from A's bow shock injected into B's open 
field line region and Af2„ a is the unknown solid angle of A's 
wind. 

The X-ray luminosity predicted by the model is in good 
agreement with our results if we assume that PSR A's wind is 
anisotropic (i.e. Afi^.A ^ 1 sr), as heuristically suggested by 
the variable illumination inducing radio and X-ray emission 
of PSRB. Within the above anisotropy assumption, the re- 
quired efficiency 77 is (poorly) constrained to <50% to match 
the observed PSR B luminosity. 

However, according to this model, PSR B's surface would 
be heated to a temperature (^500 eV) higher than the ob- 
served values of BBc and BBh. This discrepancy could be 
due to the fact that most of the predicted luminosity (BBc) 
comes from a region larger than the polar cap and consistent 
with the whole neutron star's surface. The hotter component 
is consistent in size with the polar cap and could be responsi- 
ble for the observed pulsations. 

In radio, the shape and the intensity of the pulse profile of 
PSRB vary with orbital longitude, with two bright phases 
centered around longitudes (with respect to the ascending 
node) of ^210° (bright phase 1 [bpl], orbital phase 0.11) and 



280° (bright phase 2 [bp2] , orbital phase 0. 3 1 ) at the epoch of 
first observations (Lyne et al. 2004). Because of the geode- 
tic precession of pulsars' spins and periastron advance, these 
bright phases are shifting to greater longitudes at a rate of 
a few deg yr"^ (Burgay et al. 2005). At the epoch of our 
XMM-Newton observation, the centroid of the radio bright 
phase bp2 is expected at ^300° of longitude (orbital phase 
0.36). The brightening of PSRB in X-rays appears in the 
orbital phase interval adjacent to (following) bp2 and peaks 
around the ascending node. When PSR A's wind intercepts 
PSR B's magnetosphere, it powers radio emission and it starts 
heating B's surface. The X-ray emission from B lasts for 
30± 10 min each orbit (detection significance fades outside 
this range) constraining the thermal inertial time of the neu- 
tron star in good agreement with theoretical expectations for 
external heating of a pulsar surface (Eichler & Cheng 1989). 

5.2. Constraints on X-ray emission from a bow shock 
between PSR A's wind and PSR B's magnetosphere 

The short occultation of A by B seen in the radio band 
(Lyne et al. 2004; McLaughlin et al. 2004c,b) could imply 
the presence of a dense magnetosheath enfolding the magne- 
tosphere of PSRB, although all the absorption could occur 
from within the PSR B's magnetosphere alone (Lyutikov & 
Thompson 2005). In the magnetosheath model, the relativis- 
tic wind from A collides with B's magnetic field at an equiUb- 
rium distance from B comparable or smaller than the radius 
of its light cylinder (Arons et al. 2005; Lyutikov 2004). Thus, 
in contrast with termination shocks seen in pulsar wind neb- 
ulae (PWNe) typically at distances >10^ light cylinder radii, 
a termination shock between PSR A and PSRB could probe 
the properties of a pulsar's relativistic wind at a smaller dis- 
tance from the central engine than ever studied before. PWNe 
are usually strongly radiative at high energies with typical ef- 
ficiencies Lx/EpsR = 10"^ to 10"^. Therefore, some trace 
of a peculiar PWN originating from PSR B's magnetic pres- 
sure confining the wind outflow from A should be in principle 
found as a non thermal component in our data. 

An important issue concerning such a bow 
shock/magnetosheath model is that one might expect a 
flux modulation as a function of the orbital phase owing 
to the changing view of the shock front. In particular, the 
shocked wind is expected to flow away from the head of 
the bow shock in a direction roughly parallel to the shock 
(Lyutikov 2004; Granot & Meszaros 2004). This might imply 
relativistic beaming of the radiation emitted by the shocked 
plasma, resulting in a <50% modulation of the observed 
emission as a function of the orbital phase with possible 
peaks at 90° from conjunctions. 

The upper Umits we obtained on the flux variations as a 
function of the orbital phase allow us to exclude any signif- 
icant orbital modulation ^15% in X-rays (§ 3.1). Assuming 
that most of the unpulsed flux of the system (^10%-50% of 
the total flux depending on the selected energy range) could 
be ascribed to the bow shock, flux variations should have been 
detected. In fact, the range of allowed spectral slopes of non 
thermal components in our spectral fits is not compatible with 
the "canonical" shock value F 2, although a soft post-shock 
spectrum could be explained with the presence of an unusual 
low-energy relativistic electron population, as required in the 
magnetosheath model to provide A and B's eclipses by syn- 
chrotron absorption (Arons et al. 2005; Lyutikov 2004). Nev- 
ertheless, the off-pulse spectrum of PSR A is not even com- 
patible with a simple non thermal emission, and PSRB would 
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fill most of it at least for the part of the orbit in which it is 
detected. Therefore there is not much room left for a shock 
emission component in our data. An upper limit (90% confi- 
dence level) of ~10^^ erg s~' (0.2-3 keV) corresponding to an 
efficiency of L^°^^/Ea — 2x 10~^ can be obtained by evalu- 
ating the maximum allowed luminosity of an unpulsed power 
law component with photon index ~2 in the off-pulse spec- 
trum of PSR A. 

Such a low output in shock emission compared to typical 
PWN's lJ°*/£'psr is not surprising looking at the small solid 
angle over which the wind from A is intercepted by B and 
considering that the termination shock of the wind is much 
closer than in known nebulae and that it could show a different 
phenomenology. 

Assuming that the wind energy is radiated isotropically 
from PSR A and that it is intercepted by a sphere cen- 
tered on PSRB with radius equal to its light cylinder /?ic,b, 
the resulting shock emission efficiency is L*°'=''/£psr = 
kRi^-g/4d^f{aM), where d is pulsars' separation, /(ctm) is 
the fraction of power intercepted by the shock feeding into 
the accelerated electrons (it is a function of the unknown 
magnetization parameter ctm), and k is a normalization fac- 
tor to fit observed "standard" PWNe X-ray efficiencies, where 
um 1. Kennel & Coroniti (1984a,b) provide an extensive 
discussion of the shock power function /{(Jm) that is propor- 
tional to 1/((Tm)'^^ for (Tm ^ 1 (high-(TM shocks have poor 
efficiency) and is <1 for ctm <C 1. Figure 8 shows the re- 
sulting L^^^/^PSR curve as a function of the magnetization 
parameter ctm accounting for uncertainties (hatched area) as- 
sociated with the actual size of the region intercepting the 
shock (it could be a factor 5 smaller than R\c.b according to 
pressure balance) and errors on k due to the observed span in 
PWNe efficiencies. For low values of the magnetization pa- 
rameter, the shock could theoretically account for a significant 
fraction of the X-ray luminosity, but this is not the case; the 
observed upper limit (dot-dashed line) constrains the magne- 
tization parameter to values >1 and, thus, much higher than 
those estimated for Crab-like PWNe. Such an high magneti- 
zation parameter is apparentiy inconsistent with the value of 
~0.03 predicted by the magnetosheath model (Arons et al. 
2005), and it seems to confirm most of the modern wind mod- 
els (Contopoulos & Kazanas 2002; Lyubarsky & Kirk 2001) 
predicting (Tm S> 1 . 

Granot & Meszaros (2004) proposed a model in which the 
emission from the interaction of the two pulsars is lower than 
that expected from the interaction of pulsar A's wind alone 
with the interstellar medium. Thus, PSR A could provide a 
"standard" PWN (ctm <C 1) at distances of ~10^^ cm from the 
pulsar surface. In this model, Uttle or no modulation with the 
orbital periods Pa and Pr is expected, but even in this case it 
is difficult to account for the steepness of the unpulsed spec- 
trum. Furthermore, the predicted luminosity of this model 
in the energy range (0.2-10 keV) is ~7 x 10^ erg s"\ much 
higher than our upper limit on the shock component discussed 
above. 

6. CONCLUSIONS 

The analysis of the ^5600 source photons obtained from 
our ~230 ks long XMM-Newton observation confirmed the 
pulsed detection in X-ray of the PSR A emission providing 
80% of the total absorbed flux from the system. For the first 
time. X-ray emission from PSR B is also detected with good 
confidence (^200-300 pulsed counts) around the ascending 
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Fig. 8. — Predicted X-ray emission eiBciency of a bow sliock between 
PSRA's wind and PSRB's magnetosphere as a function of tlie wind magne- 
tization parameter ctm (hatched area; see text). Tlie observed upper limit on 
the flux of a shock signature emission in our data (dot-dashed line) constrains 
Cm to values > 1 (dotted line). 

node of the orbit. 

A three-component spectral model can satisfactorily fit 
phase-resolved spectra and the energy dependence of the 
PSR A pulsed fraction. The best- fit model is composed of 
a soft power law responsible for most of the pulsed emission 
from PSR A (accounting for only ~15% of the total unab- 
sorbed luminosity at £ > 0.1 keV), a cooler (~30 eV) black- 
body likely associated to PSR B, and a hotter and fainter ther- 
mal component originating from back falling particles heating 
polar caps of either PSR A or PSRB (7?bBh ^ 100 m). Note 
that the cooler blackbody accounts for most of the bolomet- 
ric luminosity, even if its contribution in the XMM-Newton 
energy range is small. 

PSR A shows peculiar properties with respect to other 
known recycled pulsars (e.g. the softest non thermal spec- 
trum), possibly due to the particular evolutionary history of 
a double-neutron star system. The phenomenology of PSRB 
is completely different from that of any other pulsar and the 
only viable possibility to explain its X-ray emission seems 
that it is powered by PSRA's wind heating B's surface. Our 
observation is in good agreement with theoretical predictions 
by models also explaining B's radio flares (Zhang & Loeb 

2004) , although a more complex thermal emission scenario 
(colder and non uniform temperature) should be invoked. 

There are no signs (e.g. orbital flux modulation, a signifi- 
cant non thermal component in PSR A's off -pulse) of the pres- 
ence of X-ray emission from a bow shock between PSRA's 
wind and PSRB's magnetosphere (Lyutikov 2004; Arons 
et al. 2005) invoked to explain the occultation of the radio 
emission of PSR A at the inferior conjunction of PSR B. The 
upper limit on the flux of such a shock component constrains 
the wind magnetization parameter ctm of PSR A to values > 1 , 
much higher than that predicted by the magnetosheath radio 
occultation of PSR A. The absorption causing A's occulta- 
tion is likely to occur within B's magnetosphere which retains 
enough plasma to produce an eclipse (Lyutikov & Thompson 

2005) . 

Further X-ray observations of PSR J0737-3039 witii XMM- 
Newton and Chandra with reasonable exposure times (<1 
Ms) would not significantly improve PSR B detection or bet- 
ter constrain the pulsar's spectra. At gamma-ray energies, 
the absence of the strong thermal components seen in X- 
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rays could instead allow one to reveal non thermal emission 
from the pulsars' interactions and better constrain the mag- 
netospheric emission from PSR A. The empirical relation be- 
tween gamma-ray luminosity and E for known gamma-ray 
pulsars (i.e, °= Thompson 2004; Zhang & Hard- 

ing 2000) yields L^,a ~ 10^^ ergs \ a value compatible 
with the luminosity of the possibly related EGRET source 
3EGJ0747-3412 (Hartman et al. 1999) assuming the same 
distance. The gamma-ray pulsed emission of PSR A could 
be detected by the AGILE (Tavani et al. 2006) and GLAST'^ 
missions. 

Apart from PSR J0737-3039, the only other known 
X-ray-emitting double-neutron star binary (DNSB) is 
PSRJ1537+1155 (Kargaltsev et al. 2006). The X-ray spec- 
tra of PSR J1537+1155 and PSR J0737-3039 are similar, and 
their X-ray luminosities are about the same fraction of the 
respective spin-down energies. We note that, if the other 
DNSBs also have a similar efficiency, no other known ob- 
ject of this class can be detected by the currently available X- 
ray telescopes. Unlike PSR J0737-3039, in PSR J1537+1 155 
the distribution of photon arrival times over the binary orbital 
phase shows a deficit of X-ray emission around apastron, sup- 

^ See tittp://glast.gsfc.nasa.gov/. 



porting the idea that most of the emission is caused by inter- 
action of the relativistic wind from the pulsar with its neutron 
star companion. The difference between the two DNSBs can 
be explained by the smaller eccentricity of PSR J0737-3039 
or different alignments between the equatorial plane of the 
millisecond pulsar and the orbital plane of the binary (Kar- 
galtsev et al. 2006). Despite PSRJ1537+1155 being weaker 
than PSR J0737-3039, it may offer better diagnostics to dis- 
criminate emission produced by the shocked wind from "stan- 
dard" pulsar magnetospheric and surface emission, and thus, 
it is one of the most interesting systems suitable for further 
X-ray investigations of the inner pulsar magnetospheres and 
wind parameters in compact systems. 
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from our programs and Tempo. The authors acknowledge 
the support of the Italian Space Agency (contract ASl/lNAF 
1/023/05/0). A.D.L. acknowledges an Italian Space Agency 
fellowship. 
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